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Abstract 

We study soft gluon fct-resummation and the relevance of InfraRed (IR) glu- 
ons for the energy dependence of total hadronic cross-sections. In our model, 
consistency with the Froissart bound is directly related to the ansatz that 
the IR behaviour of the QCD coupling constant follows an inverse power law. 
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1. Introduction 

In this paper, we discuss soft gluon fct-resummation in the InfraRed (IR) 
region, with the aim to connect it to the energy dependence of the total 
hadronic cross-section in the high energy region. We shall make use of a 
model for total cross-sections, which incorporates in an eikonal formu- 
lation such QCD inputs as mini-jets and soft gluon fcj-resummation. This 
model has been successfully applied both to proton and photon processes: 
our aim here is to describe its physical content, without explicit reference 
to data fitting, and explain how the model incorporates a taming effect on 
the rapidly rising QCD cross-sections, thus inducing a more temperate rise 
of the total cross-section. 
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The energy behaviour of the total hadronic cross-section has been the 
focus of both theoretical and experimental enquiries for a long time. Fits 
inspired by theoretical arguments have been the subject of many debates. A 
traditional Regge-Pomeron type fit 0] such as 

(Ttotal = XS-'^ + YS% (1) 

with ?7, e > 0, presents the difficulty of not agreeing with the Froissart- 
Lukaszuk-Martin bound ji], 0, 0] 

TT 

(y total < ^ln^(s/so). (2) 
mi 

TT 

However, with ^/sq ~ 0{1 GeV), the large constant factor at the r.h.s of 
Eq. ([2]) makes possible for Eq. ([1]) to be valid only in the present region, and 
the observed behaviour to be not yet asymptotic. Still, phenomenological 
reasons run against the validity of the Regge-Pomeron type expression with 
a universal value for e, since LEP data on photon-photon indicate that the 
power, with which crJ^oJ„; rises, differs from that in pp/pp 

Currently, many models focus on QCD perturbative processes to drive 
the rise with energy, with QCD inspired models on the one hand 0, @] and 
approaches based on Reggeon calculus j^, on the other. However, it 
is still being debated how to implement this dynamics and simultaneouly 
describe both the early rise, which starts for y/s between 10 ^ 20 GeV and 



50 -7- 60 GeV, and the subsequent levelling off at higher energies . 

In this paper we shall show how both of the above features are present 
in our model for the total cross-section [l| and relate the rate, at which the 
cross-section asymptotically rises, to the IR limit of soft fct-resummation, 
thus linking directly the rise of the total cross-section to the IR region of 
QCD. We start by recalling in Sect. [2] some features of resummation and 
present, in Sect. |3l our proposal for handling the ultra soft gluon emissions 
which affect scattering at the very large impact parameter values, relevant 
to total cross-section calculations. After a brief discussion of our model in 
Sect, m we show in Sect. [5] the connection between the Froissart bound and 
our proposed eikonalized mini-jet model with ultra soft gluons. Finally, in 
Sect. |6]we discuss both the various energy scales and the constants involved 
in the model. We find that resummation of soft gluons emitted with very 
small transverse momentum kt introduces a new energy scale, which modifies 
the constant in front of the asymptotic limit derived in our model for the total 
cross-section. 
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2. Resummation and the IR limit 



The high energy behaviour of the total cross-section depends on the prop- 
erties of the scattering amplitude at large values of the impact parameter b 
in the plane perpendicular to the scattering. Recent attempts to study this 
large-6 behaviour have focused on relating Yang-Mills theories to string the- 



ories through the AdS/CFT correspondance [12, ll3|, [ij]. Our approach to 



is 



the large 6-limit of the impact picture in the eikonal representation [15 
of a more phenomenological nature: within such a picture, we exploit soft 
gluon /ct-resummation in the IR region to describe matter distribution inside 
the hadrons as they engage in hard scattering and their parton constituents 
"see" each other. We have already supplied phenomenological evidence for 
the applicability of our model to high energy scattering [if. Here we discuss 
some features of soft gluon fcj-resummation which bear on the asymptotic 
form of the total cross-section. 

We start by recalling some properties of soft photon resummation. In 
QED, the general expression for soft photon resummation in the energy- 
momentum variable Kn can be obtained order by order in perturbation theory 



16j,[l7 



as 

d^P{K) = d^K j _0_e^^— '^C-'^) (3) 

where d'^P{K) is the probability for an overall 4-momentum escaping 
detection, 

h{x,E)= d^n{k)[l - e-'^-""] (4) 
Jo 

with d^n{k) being the single soft photon differential spectrum, and E the 
maximum energy allowed for single photon emission. 

Eq. ([3]) leads to the well known, power-like, form of the energy distribution 



18l |. This is not possible for the momentum distribution, but it is also not 
necessary, since the first order expression in aqEo is adequate. On the other 
hand, for strong interactions we require the resummed, transverse momentum 
distribution, namely 

d^P{K^) = d^^±J^ j d'h e-*Kx-b-h(6,£;) 

with 

h{b,E) = [ d'^n{k)[l - e'^^-^] (6) 
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For large transverse momentum values, by neglecting the second term and 
using a constant cut-off as lower limit of integration, the above expression 
coincides with the Sudakov form factor 



IS 



Eq. p has been applied in QCD [201, M, 122, 123|, 12| with 

but its use is complicated by our ignorance of the IR behaviour of the theory. 
To overcome the difficulty arising from the IR region, the function h{b,E), 
which describes the relative transverse momentum distribution induced by 
soft gluon emission from a pair of, initially collinear, colliding partons at LO, 
is split into 

h{b,E) = coifi,b,E) + Ah{b,E), (8) 

where 

Ah{b,E) = - -^^[l- J^{bk,)]-^ In— . 9) 
S 7T kt kt 

Since the integral in Ah{b, E) now extends down to a scale /i 7^ 0, for 
/i > Aqcd one can use the asymptotic freedom expression for as{k^). Fur- 
thermore, having excluded the zero momentum region from the integration, 
Joipkt) is assumed to oscillate to zero and neglected. The integrand in Eq. ([9]) 
is now independent of b and the integral can be performed. In the range 
1/E <b < 1/A, the effective heff{b,E) is obtained by setting fi = 1/6 |2l|. 
This choice of the scale introduces a cut-off in impact parameter space which 
is stronger than any power, since the radiation function, for Nf = 4, is now 
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^eff{b,E) ^ r ln(l/&^A^) (16/25) MEVA^) 

4n(EVA2)^ ^ ' 

The remaining 6-dependent term, namely exp[— co(/i, b, E)], is dropped, a rea- 
sonable approximation if one assumes that there is no physical singularity in 
the range of integration < kt < 1/b. This contribution however reappears 
as an energy independent smearing function which reproduces phenomeno- 
logically the effects of an intrinsic transverse momentum of partons. For most 
applications, this may be a good approximation. However, when the integra- 
tion in impact parameter space extends to very large-6 values, as is the case 
for the calculation of total cross-sections, the IR region may be important 
and the possibility of a physical singularity for as in the IR region becomes 
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relevant. It is this possibility, which we exploit in studying scattering in the 
very large impact parameter region, b oo. 

3. A proposal for the IR limit in the soft gluon integral 

In this Section we discuss a phenomenological expression for the coupling 
of ultra soft gluons to the emitting quarks, and compare the resulting large-6 
behaviour with the one discussed in the previous section. 

Our choice for the IR behaviour of ^^(Q^) used in obtaining a quantitative 
description of the distribution in Eg. (E l), is inspired by the Richardson 
potential for quarkonium bound states [25], as we have proposed in a number 



of related applications 26[]. Assume a confining potential (in momentum 



space) given by the one gluon exchange term 

V{Q) = K{^^), (11) 

where K is a constant, calculable from the asymptotic form of ^^(Q^)- Let 
us choose for << the simple form 

-M') = Jg^,. (12) 

(with B a constant), so that V{Q) for small Q goes as 

V{Q)-^Q~^^'+'^\ (13) 



For the potential, in coordinate space, V{r) = J d^Q/{2TxYe'-^-''V{Q), Eq. f|T3|) 
implies 

V{r) {l/rf ■ r(2+2p) ^ C r^^P"^\ (14) 

for large r (C is another constant). A simple check is that for p equal to 
zero, the usual Coulomb potential is regained. Notice that for a potential 
rising with r, one needs p > 1/2. Thus, for 1/2 < p < 1, this corresponds 
to a confining potential rising less than linearly with the interquark distance 
r, while a value of p = 1 coincides with the IR limit of the Richardson's 
potential and is also found in a number of applications to potential estimates 



of quarkonium properties [27 . 
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Then, again following Richardson's argument, we connect our IR limit 
for UsiQ"^) to the asymptotic freedom region using the phenomenological 
expression: 

_ 127r p 



(33-2iV^)ln[l+p(^)p] 



(15) 



which coincides with the usual one-loop formula for values of kt » A, while 
going to a singular limit for small kt, and generalizes Richardson's ansatz 
to values of p < 1. The range p < 1 has an important advantage, i.e., it 
allows the integration in Eq.(l6]) to converge for all values of kt = \k±\. Using 
Eq. ( |T5i) . one can study the behaviour of h{b,E) for the very large-6 values 
which enter the total cross-section calculation and recover the perturbative 
calculation as well. The behaviour of h{b, E) in various regions in 6-space was 
discussed in 28|, both for a singular and a frozen namely one whose IR 
limit is a constant. There we saw that, for the singular case, the following 
is a good analytical approximation in the very large-6 region: 



h{b, M, 



2cf 

TT 

2cf 

IT 







dk , 2M 

In — + 2bA^P 



p) 



k^P 



dk ^ M - 



dklnf 



2 ln(2M6) + 



1 — p 



2pM; 



2p 



-2 In 



M 

aK 



1 

+ - 
p 



- M 
+ 61n — 
A 



In ^ 



IniV, 



2 ln(M6) - 
1 + 



This approximation is valid in the region b > l/{NpA) > 1/M, with Np = 
{1/pY^'^p, Cp = 4/3 for emission from quark legs and b = 127r/(33 — 2Nf). 
The upper limit of integration, M, indicates the maximum allowed transverse 
momentum, to be determined by the kinematics of single gluon emission as in 



29(1 . The above expression exhibits the sharp cut-off at large-6 values which 
we shall exploit to study the very large energy behaviour of our model. On 
the other hand, the possibility that as becomes constant in the IR 2l|, |22|, |23 
in the same large 6-limit leads to 



1 

p 

InNr 



Inf 



+ 



(16) 



h{b,M,A) = (constant) \n{2Mb) + double logs 
namely no sharp cut-off in impact parameter b, as expected. 



(17) 
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4. The Bloch-Nordsieck(BN) model for the total cross-section 



Our model for the total cross-section [Jj is a modified mini-jet model, in 
which the rise with energy is driven by perturbative parton-parton scatter- 



ing [30], tempered by an energy dependent acollinearity effect. This effect is 
due to fc^-resummation of soft gluon emission from the initial state, hereafter 
referred to as soft gluon kt- emission. The emphasis on resummation, first in- 



troduced for electron scattering by Bloch and Nordsieck [31[ , gives the model 
its name. The model is built through the eikonal representation in impact 
parameter space, so as to satisfy unitarity, and allows to implement multiple 
parton scattering and to restore a finite size of the interaction through the 
impact parameter distribution in the scattering hadrons. The details of the 
model can be found in [l|, [isl, 32], here we shall recall some aspects relevant 
to its asymptotic energy behaviour. 

In hadron-hadron scattering at a cm. energy y/s, unitarity allows to 
write a simple model for the total cross-section, namely 

atotai = 2 j d%[l - e-^'"'^(^'^)cos3?ex(&, s)] ^ 2 j (fh[l - e-"^^'^)/^] ^^g^ 

where the approximation on the r.h.s is obtained by neglecting the real part 
of the eikonal function (at the hadronic level, an acceptable approximation 
in the high energy limit) and 2Xmx{b, s) = n{b,s). The latter follows from 
a semiclassical argument relating aineiasUc to a sum of a Poisson distributed 
independent, single and multiple collisions. 

We use perturbative QCD to calculate the cross-section in order to obtain 
the average number of inelastic collisions. While implementing the QCD 
calculation, albeit approximate, we distinguish between the average number 
of collisions receiving contributions from hard physics processes and those 
from non perturbative ones, and write n(6, s) in the form 

n{h, s) = nNp{b, s) + Uhardib, s) (19) 

where the non perturbative (NP) term parametrizes the contribution of all 
those processes for which initial partons scatter with pt < ptmin, with ptmin 
a suitable low energy cut-off for the QCD parton-parton cross-section. We 
parametrize nNp{b,s), which establishes the overall normalization, and focus 
our attention on the hard term, which is responsible for the high-energy rise 
and which we expect to dominate in the extremely high energy limit. We 
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approximate this term as 



rihardip, s) = A{b, s)ajet{s). (20) 

The QCD jet cross section drives the rise due to the increase with energy of 
the number of partonic colhsions. It is calculated from the usual perturbative 
QCD expression, with DGLAP evoluted parton densities and perturbative 
partonic differential cross-sections. In Eq. fl20|) . A{b,s) is the overlap func- 
tion which depends on the (energy dependent) spatial distribution of partons 
inside the colliding hadrons, averaged over the densities [H, [s^] ■ Before dis- 
cussing this function we shall examine the energy behaviour of the mini-jet 
cross-sections. 

In the ^/s >> Ptmin limit, the major contribution to the mini-jet cross- 
sections comes from collisions of gluons carrying small momentum fractions 
Xi^2 << I5 a region where the relevant PDFs behave approximately like 
powers of the momentum fraction with J ~ 1.3 [33*1. This leads to the 
asymptotic high-energy expression for ajet 



1 - 



Pt mill 



t mm J 



(21) 



where the dominant term is a power of s. Fits to the mini-jet cross-sections. 



obtained with different PDF sets 34 confirm the value e = J — 1 ~ 0.3 



Such energy behaviour as in Eq. (1211) is at odds with the gentle rise 
of the total pp and pp cross-sections at very high energy, described rather 



as Ins, In s [11[ or s power [2]. However, we shall see that a proper 
implementation of other QCD processes can modify this strong rise. To do 
so we now examine the energy behaviour of the next component of our BN 
model, the impact parameter distribution. 

We have identified soft gluon /cj-emission from the colliding partons as 
the physical effect responsible for the attenuation of the rise of the total 
cross section. These soft emissions break collinearity between the colliding 
partons, diminishing the efficiency of the scattering process. Their number 
increases with energy and thus their contribution remains important, even at 
very high energy, influencing matter distribution inside the hadrons, hence 
changing the overlap function, which is proposed to be the Fourier transform 
of the previous expression for the soft gluon transverse momentum resummed 
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distribution, i.e,. we put 



(22) 

The integral in qmax) is performed up to a value gmax, which is linked to 
the maximum transverse momentum allowed by kinematics of single gluon 
emission 29|]. In principle, this parameter and the overlap function should 
be calculated for each partonic sub-process, but in the partial factorization of 
Eq. fl20l) we use an average value of qmax obtained by considering all the sub- 



processes that can happen for a given energy of the main hadronic process [28 
The energy parameter q^ax is of the order of magnitude of ptmin- For present 
low— X behaviour of the PDFs, in the high energy limit, qmax is a slowly vary- 
ing function of s, starting as In s, with a limiting behaviour which depends 



on the densities 35|. From Eqs. flTBl) and fl^T]) one can estimate the very 
large s-limit 

nhard{h,s) = ABN{b,s)a,et{s,Ptmin) ~ ^(s)^^^^'^— Vi (-)^ (23) 

So 



and, from this, using the very large 6-limit, 



Uhardib, s) ~ Ao{s)a,e-^''^^''{-r (24) 

So 

with Ao{s) oc and with a logarithmic dependence on qmax, i-e. a very 
slowly varying function of s. We also have 

A ^ A(6, .) = A{ [H2qmaxis)b) + -^]V/'^ (25) 

47r(l — p) 1 — p 

In the next section, we shall see how the two critical exponents of our model, 
namely the power e with which the mini-jet cross-section increases with en- 
ergy and the parameter p dictating the IR behaviour of the QCD coupling 
constant, combine to obtain a rise of the total cross-section in agreement 
with the In^ s limitation imposed by the Froissart bound. 



5. Ultra soft gluons in the IR limit and the asymptotic limit of 
the total cross-section: the Froissart bound 

In this section we consider the very large s-limit of the total cross-section, 
in the approximation that all constant (or decreasing ) terms in the eikonal 
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function can be neglected, thus studying only the QCD effects from mini-jets 
and soft /^^-resummation. We find a link between the infrared behaviour of 
the ultra soft gluons and the asymptotic Froissart-like behaviour of the total 
cross-section and discuss it. 

Let us consider the total cross-section in the eikonal representation at 
very large asymptotic energies. At such large energies that n^p << Uhard, 
the total cross-section in our model jl[ reads 

oo 

2ri _ „-nhard{b,s) /2-I 



Jo 

We consider the asymptotic expression for ajet at high energies, which grows 
like a power of s, and ABN{b,s), which was obtained through soft gluon 
resummation, and which decreases in 6-space at least like an exponential 
(1 < 2p < 2). In such large-6, large- s limit, we can write 

nw = 2C(s)e-(^^)''' (27) 
where 2C(s) = Ao{s)ai{s/soY. The resulting expression for ar is 

/•oo _ „ 

27r/ rf62[l-e-^W'=-'"''] (28) 
Jo 



axis 



With the variable transformation u = [Ab)"^'^, and neglecting the logarithmic 
6-dependence in A by putting b = 1/A, Eq. fl28l) becomes 

aj,{s) ^ — i- r duu'/P-'[l - e-^(^)^-"] (29) 
P A^ Jo 

Notice that, as s oo, C{s) also grows indefinitely as a power law. This 
means that the quantity between square brackets I{u, s) = 1 — e~'-^('*)^ " has 
the limits I{u, s) — > 1 at m = and I{u, s) — > as u = oo. Calling uq the 
value at which /(mq, s) = 1/2 we then put I{u, s) ~ 1 and integrate only up 
to Uq. Thus 

AVt(s) ^ (— ) / " duu^ = (30) 
P Jo 



and since, by construction 



Uo = H^]^e\ns (31) 
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we finally obtain 

ar^%[e\n^Y'^ (32) 

to leading terms in Ins. We therefore derive the asymptotic energy depen- 
dence 

ar ^ [eHs)f/^'^ (33) 

apart from a possible very slow s-dependence from A^. The same result is 
also obtained using the saddle point method. 

This indicates that the Froissart bound is saturated if p = 1/2. We shall 
now show that, in our model, analyticity demands p > 1/2 and thus that, 
no matter how fast the mini-jet cross-section may grow with energy, the 
Proissart bound is always satisfied. 

The requirement that p > 1/2 follows from analyticity arguments in the 
complex ^s-plane, where Zg = cosOg = 1 + 2t/(s — 4m^) ^ 1 + 2t/s is the 
cosine of the s-channcl scattering angle for the equal mass case. Basically, this 
limitation comes from asking that for large values of the impact parameter 
space b, the eikonal be such to decrease at least like an exponential, i.e. to 
have at least 

Imx{b, s) ^ e-''^ (34) 

where to > is the boundary of the Lehmann ellipse on the real axis in the 
Zg plane. To see the argument, consider the elastic scattering amplitude for 
a process a + b ^ a + b, normalized in such a way that 

atot = 2nImF{s, t = 0) (35) 

with 

F(s, t)^ij rf'b[l - e^>^(^'^)] Jo(^'\^) (36) 

While in the s-channel physical region, and for equal mass particles, s > 
4m^ and t < 0, for the Lehmann eUipse F(s, t) is analytic for t < to < jJ^^ 
where // is some hadronic mass (e.g. twice the pion mass, as the smallest 
mass exchanged in the t-channel). The actual value is unimportant, what is 
needed is to > 0. Then, for < i < io the argument of Jo becomes imaginary 
so that 

F(s, t>Q)=i J d^h[l - e'^^''''%{bVi) (37) 

where lo is the Bcssel function of the second kind. The asymptotic behaviour 
of the Bessel functions is such that Io{y) for y real, grows exponentially as y 
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becomes large. The integral at the r.h.s. of the above equation has to exist 
up to values of t = to and for fixed (0 < t < to)- Since, for large 6- values, 

loibVi) ~ e'^, (38) 

for the integrand to be finite the imaginary part of the eikonal function, 
Tmx{b, s) must go to zero at least like an exponential, i.e. at least 

Imxib, s) e-^"^ (39) 

Now, let us return to the expression for Imxib, s) from the BN model, where 
Imx{b, s) ~ e"'^*^*''^^ up to exponential accuracy. Using the large-6 behaviour 
of the function h{b,s), derived in previous papers and reproduced in the 
previous section, we see that 

Imxib, s) ~ e-'^^'''") ~ e"^^^'" (40) 
whence follows that p > 1/2. 



6. About scales and parameters in the BN model 

The model we have described contains different scales. As mentioned 
earlier, in this paper we only study the rising part of the cross- sect ion, whose 
behaviour is proposed to come from processes for which the outgoing partons 
have transverse momentum pt > ptmin — 0{1 GeV). In our model, ptmin is 
the scale which separates perturbative scattering processes from everything 
else. Soft gluon emission introduces two more scales, namely A and qmax- 
The latter is the maximum transverse momentum in the integral for soft 
gluon emission, it is of order ptmin and plays the role of the energy scale E 
which appeared in QED radiative correction factors. Thus soft gluons satisfy 
the condition 

kt < Qmax ^ Oiptmin), (41) 

since most of the parton-parton cross-section is peaked at pt = Ptmin- The 
next scale A ~ 0{Aqcd) separates the region of ultra soft gluons from the 
rest. This region was originally neglected, on the basis that gluons with 



\k±\ <C A would see the hadron as a point-like object [36|] and such emissions 
would have a small probability, because of colour screening. This argument 
is appealing, and similar to the one mentioned in Sect. [21 but in our opin- 
ion, there is no compelling theoretical reason to assume that ultrasoft gluon 
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emission in high energy reactions has low probabihty. This argument could 
be applied to an isolated hadron, but not to high energy hadronic scatter- 
ing described through the scattering of partons, where soft gluon emission is 
stimulated by QCD interactions. It is through this interaction that we can 
expect the transition between hadrons and quarks to arise. A singularity in 
the infrared region would indeed provide a cut-off to separate quarks from 
hadrons and lead to such transition. This is the rationale behind going into 
the zero momentum region, and enter it with a singular confining coupling 
between ultra-soft gluons and the quark current. In our model, these ultra 
soft gluons are important for the extremely large impact parameter values, 
which enter eikonal formulations of the total cross-section at very large en- 
ergy. For processes where such large-6 values do not play a role, this region 
may be irrelevant, though. 

The ultra soft gluon distribution which we have introduced depends on the 
parameter p which regulates the infrared region of the soft gluon integral of 
Eq. ([7]). In we have used the value p = 3/4 but this is a phenomenological 
value which was obtained after performing various averages over the PDFs. 
A determination of the actual value of p is beyond the scope of this paper, 
except for the fact that p < 1 for the soft gluon integral to converge and 
p> 1/2 for analyticity of the scattering amplitude. 

One can evaluate the coefficient of the (In s)^/^ term in Eq. (133|) as given 

by 

^ 271 , 27(1 -pfe 

mr A ^ M[l + (l-p)ln(2g^,,./A)]^ ^ ^ 

for Nf = 3, as appropriate for the total cross-section limit. In the approxi- 
mation in which the Qmax term is neglected, Eq. (l42l) gives ctt ~ vr/m^ln^s 
for p = 1/2, A = 100 MeV and e = 0.3. However, Qmax, which provides 
an extra dynamical scale, cannot, in general, be neglected: using Eq. 0421) 
with Qmax — 1 GcV , the above equation gives C ~ O(0.1)7r/m^, namely a 
constant which is one order of magnitude smaller than in the case of the 
actual Froissart bound j^. 

We note that in the eikonal model as used here, one can derive the limit 
for the inelastic cross-section 37| following the same steps as above obtaining 
a constant reduced by a factor 2, as in a black disk model. 
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Conclusions 

Using an eikonal mini-jet model for the total cross-section we have shown 
how soft gluon fcj-resummation in the IR region can reduce the strong power- 
like rise due to the minijet cross-section. We have found that this model will 
always satisfy the bound dtotai < In^ s provided that the infrared behaviour 
of as reflect a rising one gluon exchange potential: for a potential rising 
like r^P-i the total cross-section is limited by an asymptotic behaviour ~ 
(In sY^P, with 1/2 < p < 1. This establishes the connection, in our BN 
model, between confinement and the satisfaction of limitations imposed by 
the Froissart bound. 
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